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ABSTRACT 

We present an analysis of a high resolution (R ~ 50, 000) optical spectrum of 
the central region of the proto-planetary nebula CRL 2688. This object is thought to 
have recently moved off the AGB, and display abundance patterns of CNO and heavy 
elements that can provide us with important clues to understand the nucleosynthesis, 
dredge-up and mixing experienced by the envelope of the central star during its AGB 
stage of evolution. The analysis of the molecular features, presumably originated from 
the circumstellar matter provides further constraints on the chemistry and velocity of 
the expanding shell, expelled as a consequence of the strong mass loss experienced by 
the central star. 

We confirm that the central star shows a spectrum typical of an F-type super- 
giant with T cS = 7250 ± 400 K, log g = 0.5 and [Fe/H]= -0.3 ± 0.1 dex. We find 
that the abundance pattern of this object is characterized by enhancements of Carbon 
([C/Fe]= 0.6 ± 0.1), Nitrogen ([N/Fe]= 1.0 ± 0.3) and Na ([Na/Fe]= 0.7 ± 0.1), simi- 
lar to other previously known carbon-rich post-AGB stars. Yttrium is also enhanced 
while the [Ba/Y] ratio is very low (—1.0), indicating that only the light s-process ele- 
ments are enhanced. The Zinc abundance is found to be normal, [Zn/Fc]= 0.0 ± 0.3, 
suggesting that there is no depletion of refractory elements. The Ha, Nai and Ki 
resonance lines show prominent emission components, whose helio-centric radial ve- 
locities are offsetted by —41 ± 3 km s _1 relative to the photospheric metal-absorption 
lines. The molecular features of C2 and CN also show emission components, whose 
velocities are consistent with the emission components of the Ha, Nai and Ki lines. 
On the other hand, their absorption components are more highly blue shifted than 
the corresponding emission components, which suggests that the regions where the 
emission and absorption components arise are expanding at different velocities. 

Key words: stars: AGB and post-AGB - stars: evolution -stars: abundances- stars: 
circumstellar matter -stars: individual (CRL2688). 



1 INTRODUCTION 

CRL 2688 (V1610 Cyg, Egg Nebula) is a bipolar 
proto-planetary nebula (PPN). The central star is 
thought to be a F5Iae post-AGB supergiant. Initial 
studies of CRL 2688 were mostly based on low res- 
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olution spectrophotome tric and spec tr opolarimetric 



data and were made by Nev et al. Jl975 ). iForrest et al.l 
ll 19751). ICrampton. Cowley fe Humphreys! (|l975l l and 
ICohen fc Kuhil (|l9T7l . Il980h . The low resolution spec- 
trophotometric observations revealed emission from C2, 
C3 and SiC2- The Swan bands o f C2 are also s een i n 
emission from the bipolar lobes (|Cohen fc Kuhil Il980h . 
Because the object is faint (V = 13.5) with significant 
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circumstellar reddening its high resolution spectroscopy 
was not carried out until the year 2000. The first high 
res olution (R = 1 5 ,000) optical spectra were analyzed 
by iKlochkova et "all boOOh . They found that the central 
star is a F5Iae post-AGB supergi ant with overabundance 
of carbon and s-process elements. IKlochkova et al.l (2000) 
derived T c s = 6500 K, \ogg = 0.0, microturbulent velocity 
= 6.0 km s^ 1 and [Fe/H]= -0.59. 

High-spatial resolution Hubble Space Telescope imaging 
revealed the detailed morphology of this object; a flattened 
dust cocoon obscuring the central star as well as extended 
bipolar nebulae exhibiting numerous concentric arcs crossed 
by a pair of bea m-like structures a t the both sides of the 
equatorial plane l|Sahai et al.l 1 19981 ). These arcs are sug- 
gested to be associa ted with repeated mass-loss events d ur- 
ing the AGB phase i|Sahai et al.lll998l ; iBalick et af1l2012l ). 

Recently, IWesson et all (|20ld ) studied the Herschel 
Space Observatory - SPIRE FTS spectra of CRL 2688 cov- 
ering the wavelength range 195 to 670 /jm. They found the 
far infrared spectrum of CRL 2688 to be very complex with 
about 18 different species of complex molecules. They also 
detected water in CRL 2688 for the first time. 

CRL 2688 shows small amplitude light vari ations with 
a period of about 90 days l|Hrivnak et al.l 120101 ) and some 
of the absorption and emission lines seems to show a radial 
velocity variability with time (|Klochkova et all [2000) . The 
variability indicates that the central star as well as the cir- 
cumstellar material is undergoing rapid variations in their 
dynamical and chemical properties possibly due to pulsation 
and mass- loss. 

Mechanisms that give rise to observed complex mor- 
phology, chemistry and variability of CRL 2688 as well as 
other post-AGB stars are still far from being completely 
understood. Identifying constraints on these mechanisms is 
important for understanding the chemical evolution of our 
Galaxy since low-to-intermediate mass stars on the AGB 
phase of evolution are thought to be mainly responsible for 
synthesizing and contributing with CNO and s-process ele- 
ments to the interstellar medium via mass loss. 

In this paper we present an analysis of a much higher 
resolution spectrum of CRL 2688. 



2 OBSERVATION AND REDUCTION 

A high resolution (R ~ 50, 000) spectrum of CRL 2688 was 
obtained with the Utrecht Echelle Spectrograph (UES) at 
the 4.2 m William Herschel Telescope on July 14, 2001. The 
slit of 1 arcsec widths was positioned at the central star 
region. The wavelength range from 4300 to 9000 A was ob- 
served with an exposure time of 1200 sec. The data were 
reduced with standard IRAF routines. The signal to noise 
ratio of the reduced spectrum per resolution element varies 
from 30 to 70 in this wavelength range. 

Equivalent widths (EW) of atomic absorption lines were 
measured by fitting Gaussian to each feature using a linelist 
described below. 



3 DESCRIPTION OF THE SPECTRUM 

The reduced spectrum shows various photospheric absorp- 
tion lines and molecular lines. A portion of the observed 
spectrum is shown in Figure[T] Absorption lin es that are con- 
sidere d in t he EW measurements t aken from lRevniers et al.l 
l|2004h and IKlochkova et all (|2000h are marked. 



3.1 Atomic lines 

The obtained optical spectrum of CRL 2688 exhibits nu- 
merous atomic absorption lines of metals, most of which are 
too weak to be used in the abundance analysis, as illus- 
trated in Figure [T] Many of the relatively strong lines like 
H a and Ba 11 show a blueshifted emission component super- 
imposed on the absorption component as shown in Figure [2] 
The presence of emission component in the H a violet wing 
is in agreement with the suggestion from a spectropolari- 
metric observation of IKlochkova et all |2004). The observed 
features indicate mass loss from the central star. 



3.2 Molecular lines 

The molecular lines in this star are likely originated from its 
circumstellar envelope, because such molecules would be de- 
stroyed in its photosphere with T e g ~ 7000 K. The spectral 
regions dominated by molecular bands are shown in Figure 
E] and |H 

Strong emission bands of the C2 Swan sys- 
tem are identified consistent with pr evious studies 
ilCrampton. Cowley fc Humphreys! 1 19751 ; ICohen fc Kuhil 
1 19801 ; IKlochkova et alj|200d . 12004! ) . The particularly strong 
band at 5635. 5 A of C2(0,l) i s show n in the middle panel 
of Figure [3] ICohen fc Kuhj (|l977l ) showed that the C 2 
emission bands are prominently seen only in unpolarized 
light. The negligible polarization for the C2 Swan emission 
bands has also reported by the study of IKlochkova et al.l 
(2004) with a much higher-resolution spectropolarimetry, 
suggesting that the observed C2 emission may come directly 
from the line forming region. 

We confirm the presence of molecular absorption bands 
of the C2 Phillips system and the CN red system, as shown in 
Figure 31 The detection of these molecular absorptio n bands 
had previously been reported bv lBakker et al.l (|l997l ). Addi- 
tionally, some of the features in C 2 Phillips (2,0) and (3,0) 
bands are shown in emission, whose central wavelength is 
less blue shifted than the corresponding absorption compo- 
nent. A few transitions in the CN (2,0) band also appear to 
sho w emission component s. 

ICohen fc Kuhil (|l980T) have suggested the presence of 
the Merrill-Sanford bands of SiC2 in absorption, strongest 
at 4977 A, particularly in the region of the nebular lobe. 
Our spectrum does not show this feature in agreement with 
the apparent lack of this feat ure in the study of th e central 
region of CRL 2688 made bv lCohen fc Kuhj (|l980l ). 



3.3 Nal and Kl resonance lines 

Figure [5] shows portions of the spectrum around Nal and K I 
resonance lines. Nal lines are expected to have contribution 
from stellar photosphere, circumstellar envelope and inter- 
stellar matter. However as we will indicate in Section 13.41 
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Figure 1. A portion of the obtained spectrum of CRL2688. Positions of metal-lines used in the EW determination are indicated. 



[htbp] 



6 

o 
2 









1 1 


H a 6563 A 






Ba II 6142 A . 






Ba II 6497 A . 


1 .... 1 ■ . 




... 1 .... 1 , , 



-S00 -100 100 

Velocity km s -1 



Figure 2. Line profile for Ha and Ba II lines. Only the core 
region of the line profile is shown for the Ho line. All lines are 
shifted in velocity relative to the photospheric one. 



contribution from the interstellar absorption is likely small 
for this object. The observed Nai Dl and D2 lines show 
a complex structure showing multiple components. Each of 
the Dl and D2 line shows at least two emission components 
('el' and 'e2' in the top panel of Figure [5]). The stronger 



component ('el') is blueshifted while the weaker component 
('e2') is located close to the restframe wavelength. The K I 
lines of this star are likely of circumstellar origin since K 
would be mostly ionized in the photosphere. Similar to the 
Nai lines, the Kl line shows both absorption and emission 
components ('al','a2' and 'el'). The absorption component, 
which is at the bluer side of the emission feature, appear to 
show at least two velocity components. 

Signatures of multiple components of Nai and 
Kl absorption features have been report ed in other 
post-AGB stars, like e.g . in HD 56126 l|Bakker et all 
1996; ICrawford fc Barlowl l2000h . Based on a very high- 
resolution (R ~ 90000 0) observations of this object, 
ICrawford fc Barlowl ((2000) reported that the K I absorption 
feature of this object shows multiple components, some of 
which have a counterpart in C2 molecular features. They 
suggested that the individual components are associated 
with distinct shells of circumstellar matter around this ob- 
ject. Higher resolution spectroscopy is clearly needed to ad- 
dress the nature of the shape of Na I and K I lines in CRL 
2688 by resolving their intrinsic profiles. 

iKlochkova et al.l (120041 ) reported that the Nai emission 
component is likely unpolarized, suggesting that it is came 
directory from the observed lobe region. The presence of the 
emission components in both Nai and Kl lines as seen in 
CRL 2688 has also been observe d in a number of o ther C-rich 
(and O-rich) post-AGB stars l|Luna et alJ l^S). although 
some of the similar class of objects like IRAS 06530-0123 d o 
not seem to show these emissions (Hrivnak fc Reddvll2003l ). 
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Figu re 3. Th e obse rved spectrum of the C2 Swan system. The C2 emission features are indicated based on the identification by 
ICohen fc KuhH dl977l l. 
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Figure 5. Nal and Kl lines. Emission and absorption compo- 
nents for which radial velocities of the line center have been mea- 
sured are indicated by dotted lines. Their restframe wavelengths 
are indicated by arrows. 



3.4 Diffuse interstellar bands 

Equivalent widths of diffuse interstellar bands (DIBs) have 
been calibr ated to estimate th e line-of-sight reddening 
(E(B — V)) (|Friedman et aljfeoill ). In the present spectrum 
of CRL2688, most of the reported strong DIBs are too weak 
or likely contaminated by neighboring spectral features. The 
well calibrated DIB centered at 5780.5 A may also be con- 
taminated by photospheric absorption lines. Even though, 
it may be useful to estimate an upper limit of E(B — V) 
toward the line-of-sight of CRL 2688. Note that a Doppler 
shift of the observed DIB is not consistent with that of the 
molecular absorption and emission lines, which indicates the 
DIB is originated from interstel lar matter rather t han from 
circumstellar envelope. Indeed, iLuna et al l l|2008h did not 
detect signatures of the diffuse bands of circumstellar origin 
for a representative sample of Galactic post-AGB stars. 

The equivalent width of the fea ture near 5780.5 A is 67 
mA. According to the calibration of lFriedman et ail (|201ll ). 
this corresponds to E(B — V) = 0.1, suggesting that the in- 
terstellar reddening is very low at the line-of-sight of this ob- 
ject. Although CRL 2688 is located at low Galactic latitude 
at (I, b) = (80.166, —6.502), where integrated red dening is 
E(B - V) = 0.28 according to lSchlegel et all (|l998f l, the low 
E(B — V) value may be reasonable since t he distance to th is 

^ l2006h. 



object was estimated to be only 420 pc dUeta et al 



On the other hand, photometric observation of Nev et all 
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Figure 4. The observed spectru m of the C2 Phillips system (top) and the CN red system (bottom). Identification of molecular features 
are made based on the linelist of lBakker et al.l jl996l ). The blended features are indicated as 'blend'. 
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(1975) suggests that the extinction toward this object in 
visual wavelengths is at least 3 mag. If the interstellar red- 
dening is very low as suggested by the present analysis, most 
of the extinction toward the CRL 2688 is due to the circum- 
stellar matter and the dusty disk. 



4 RADIAL VELOCITY 

Radial velocities from unblended metal absorption lines have 
been measured by fitting Gaussian to individual features. 
For the spectral lines containing both emission and absorp- 
tion components, we have only used a core region of each fea- 
ture to determine its central wavelength with the IRAF splot 
routine. In the following, all radial velocities are those cor- 
rected to helio-centric values (Vq). Table [1] summarizes he- 
liocentric radial velocities obtained from the observed spec- 
tral features. 

The mean of radial velocities from unblended metal ab- 
sorption lines is —14.4 ± 0.6 km s" 1 , which represents the 
radial velocity of the central star. The standard deviation of 
the radial velocities from individual lines is a — 2.7 km s _1 . 
In the following, an error in a velocity measured from a single 
line is taken to be equal to this value. The Vq measured from 
our spectrum is in good agr eement with that obtained by 
iKlochkova etafl |200d . l2004h using metal- lines in the spec- 
trum, altho ugh the observed spectrum by IKlochkova et al.l 
( |2000L l2004h is that of the northern lobe of CRL 2688, while 
the present study is based on the spectrum of the central 
star. 

Emission components of Ha, Nal ('el' in Figure[5]) and 
Kl lines all show similar radial velocity with a mean —55 
km s _1 with respect to their rest frame wavelength. This 
value is smaller in magnitude by more than ~ 20 km s _1 
fro m that obtained from the emission co mponent of Hq 
bvlKlochkova et al.l (|2000T ) . As reported by IKlochkova et al.l 
(2000), the emission components also displayed some vari- 
ability in time. Therefore, the observed difference in the ra- 
dial velocity may be attributed to both spatial and time 
variations possibly due to the circumstellar outflow. The 
weaker emission components in the Nal Dl and D2 lines 
('e2' in Figure [5]) marginally agree with lab wavelengths. 
Therefore, it is likely that the e2 components in the Nal Dl 
and D2 lines are due to night sky emission and/or there may 
be some contribution from CRL 2688. 

The absorption component in the H a line shows the ra- 
dial velocity of — 10.6±2.7 km s _1 . This value is 3.8 km s _1 
smaller in magnitude than that from metal lines. This dif- 
ference could be attributed to the presence of emission line 
in its violet wing, modifying the true central wavelength of 
the absorption feature. Similarly, Nal Dl and D2 lines show 
a broad absorption feature ('a2' in Figure [5]) with a violet 
wing. The radial velocity of this feature is —10.6 and —10.0 
km s -1 for the Dl and D2 lines, which are similar to the 
absorption component of Hq. Another absorption compo- 
nents ('al' in FigureJS]) of the Nal lines have radial velocity 
of —69.8 km s" 1 . This value is close to the absorption com- 
ponent ('a2' in Figure [5]) of the Kl line. 

The emission components of individual features in C2 
Phillips (3,0) and (2,0) bands show radial velocity of — 53. 3± 
0.9 and — 51.7±1.0 km s" 1 , respectively, that are marginally 
consistent with those of the Ha, Nal and Kl lines. For the 



Table 1. Radial velocities 



Features 


Vq (km 






absorption 


emission 


Metal lines 


-14.4 ±0.6 




H(# 


-10.6 


-55.0 


NalDl 


-69.7, -10.6 


-54.8, 12.7 


Nal D2 


-69.9, -10.0 


-55.1, 12.3 


Kl 


-77.15, -69. 67^, 


-55.2 


C 2 Phillips (3,0) 


-67.2±0.2(3£" 


-53.3 ±0.9(5) 


C 2 Phillips (2,0) 


-67.5 ± 0.6(3) 


-51.7 ± 1.0(5) 


CN red (3,0) 


-65.5 ± 0.2(3) 




CN red (2,0) 


-66.8 ± 0.5(4) 





a Errors in Vq measured from a single line are 
assumed to be 2.7 km s — 1 , which is equal to the 
scatter (c) in Vq from the metal lines. 
^ Deblending and fitting were not successfully 
made. The values for approximate local minima 
are given. 

c Number of lines used to measure the mean Vq 
for each band. 

CN red system, the radial velocities for the different emission 
components are difficult to measure because the emission 
features are weak and appeared to be affected by blending. 

On the other hand, the absorption components of the C2 
Phillips bands have a mean radial velocity of —67 km s _1 , 
which is more highly blue shifted than the corresponding 
emission components, while the absorption components of 
the CN red system show a similar radial velocity of —66 km 



5 ABUNDANCE ANALYSIS 

The abundance analysis i s performed using a LTE abun- 
dance analysis code as in lAoki et all (120091) together with 
Kuru cz NEWODF model atmosphere (|Castelli fc Kuruczl 
2003) using the EWs measured for the atomic absorption 
lines. In order to exclude the features with erroneous fitting, 
we restrict the lines used in the abundance analysis to those 
with a FWHM in a range 0.4 to 1.0, where typical FWHM 
of the metal absorption lines in our spectrum is 0.7 A. We 
also exclude very strong lines with log(EW/A) ^ —4.4. The 
measured EWs with these criteria are used in the abundance 
analysis and are listed in Tabled 

5.1 Atomic data 

We ha ve a dopted a linelist com piled from IKlochkova et al.l 
( 2000) and lRevniers et al.l (|2004l ). The adopted log gf values 
and measured equivalent widths are summarized in Table [2] 

5.2 Atmospheric parameters 

The effective temperature (T e ff) of CRL 2688 was derived 
using the Fe I excitation equilibrium, in such a way that 
the trend of the Fe abundances from individual Fe I lines 
with their excitation potentials is minimized within 0.04 dex 
eV _1 , which is the typical error of the slope. This yields 
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Table 2. Equivalent widths 



Elem A logs/ X EW Rcf/^ 

(A) (dcx) (cV) (mA) 



C I 


ATTVL AA 


-2.19 


7.49 


136.29 


R 04 


C I 


0039.06 


-1.77 


7.95 


98. 71 


R04 


C I 


6013.16 


-1.16 


8.65 


94.61 


R.04 


C I 


6014.83 


-1.58 


8.64 


63.62 


R04 


C I 


7113.18 


-0.93 


8.65 


100.95 


K00 


N I 


7468.31 


-0.17 


10.34 


147.93 


R04 


O I 


6158.19 


-0.31 


10. 74 


119.67 


R.04 


Na I 


ODSZ.DO 


-0.67 


2. 10 


43. 70 


R 04 


Na I 


5688.21 


-0.37 


2. 10 


117.00 


R 04 


Si II 


5056.02 


0.31 


10.03 


172.29 


R.04 


Si II 


5978.93 


-0.06 


10.07 


81.26 


K00 


S I 




-1.00 


7.87 


21.87 


K00 


S I 


6743.58 


-0.56 


7.87 


36.05 


R.04 


S I 


6757.16 


-0.20 


7.87 


65.03 


R04 


Ca I 


5857.46 


0.24 


2.93 


72.66 


K00 


Ca I 


6439.07 


0.39 


2.53 


111.29 


R04 


Ca I 


6462.57 


0.31 


2.52 


102.07 


R04 


Ca I 


6717.69 


-0.52 


2.71 


50.41 


K00 


Sc II 


C A A A A 

564U.99 


-1.01 


1.50 


194. 12 


K00 


Sc II 


5667.15 


-1.21 


1.50 


112.34 


R.04 


Sc II 


5669.04 


-1.09 


1.50 


186. 12 


K00 


Sc II 


C O A OA 


-1.01 


1.51 


188.86 


K00 


Ti II 


lf7nO c o 

4 < y».oz 


-2.67 


1.08 


158.12 


R 04 


Ti II 


5418. 8U 


-2.17 


1.58 


178.38 


K00 


Cr II 


4812.35 


-1.80 


3.86 


153.08 


R 04 


Cr II 


5246.77 


-2.47 


3.71 


107.48 


R 04 


Cr II 


CJ/IA A^> 

4.5 


-2.62 


3.76 


86. 20 


K00 


Cr II 


eont o T 

5305.8 f 


-2.08 


3 . 8 3 


136.25 


R 04 


Cr II 


5308.43 


-1.81 


4.07 


112.35 


R04 


Cr II 


5334.87 


-1.89 


4.07 


196.78 


R.04 


Cr II 


54U ( .60 


-2.15 


3 . 8 3 


102.92 


R 04 


Cr II 


C /I OA AO 


-2.46 


3.76 


110.42 


R 04 


Mn I 


A TO O ,10 


0.04 


2.30 


66. 73 


R.04 


Fc I 


4966.09 


-0.84 


3.33 


62.76 


R04 


Fe I 


5049.82 


-1.35 


2.28 


50.84 


R04 


Fc I 


ronl Q A 

OoUJ.oU 


-0.88 


3.28 


79. 74 


K00 


Fc I 


t /l O T 

5364.8 f 


0.23 


4.45 


67.96 


R.04 


Fc I 


C O T /IT 

536 (.41 


0.44 


4.42 


84. 23 


R 04 


Fc I 


5393.17 


-0.91 


3.24 


86.78 


K00 


Fc I 


5405.77 


-1.84 


0.99 


163.20 


K00 


Fc I 


5434.52 


-2.12 


1.01 


97.56 


R04 


Fc I 


5445.04 


0.04 


4.39 


55.27 


R04 


Fc I 


5446.92 


-1.93 


0.99 


162.08 


K00 


Fc I 


6024.06 


-0.06 


4.55 


71.55 


R04 


Fc I 


6230.72 


-1.28 


2.56 


28.13 


R04 


Fc I 


6400.00 


-0.29 


3.60 


59.87 


R04 


Fc II 


5991.37 


-3.56 


3.15 


151.98 


R04 


Fc II 


6084.10 


-3.80 


3.20 


108.07 


R04 


Fc II 


6416.92 


-2.85 


3.89 


195.32 


K00 


Fc II 


6432.68 


-3.71 


2.89 


144.65 


R04 


Fc II 


7711.71 


-2.74 


3.90 


181.66 


K00 


Ni I 


5035.37 


0.29 


3.63 


87.94 


R04 


Zn I 


4810.54 


-0.17 


4.08 


28.38 


R04 


Y II 


5289.81 


-1.85 


1.03 


46.98 


R04 


Y II 


5728.89 


-1.12 


1.84 


168.35 


K00 


Ba II 


5853.67 


-1.00 


0.60 


228.35 


K00 


a Reference for 


log gf. 


R04: 


iReyniers 


ct al. 



(2004), K00: IKIochkova et al.1 §000) 



T c ff = 7250 ± 400 K. This va l ue is h igher than the previous 
estimate by IKIochkova et al.1 (|200Ch . (T e ff = 6500 K with a 
typical uncertai nty of 200 K) but in good agreement with 
that obtained bv lBakker et al.l (| 19971 ) (T efi = 6900 K) within 
the quoted error. The surface gravity (logg) is estimated 
from the Fe I and Fe II ionization balance. The model at- 
mosphere with log g — 0.5 dex, which is a lower limit of log g 
in the available model grids with T e g > 7000 K, yields rea- 
sonable agreement between the Fe abundances based on the 
neutral and ionized species. Finally, the micro-turbulence 
velocity £ is estimated by minimizing the Fe I abundance 



versus equivalent widths correlatio n. This results in £ = 6. 5 
km s _1 which agrees with that from IKIochkova et aLlfcoOOt ). 

5.3 Abundances 

Adopting the atmospheric parameters estimated above, we 
have obtained abundances of C, N, O, Na, a elements, 
Fe-peak elements and neutron capture elements from the 
measured equivalent widths. For heavy neutron-capture el- 
ements, spectral synthesis have been applied, which gives 
upper limits in the abundances. In the following, we have 
taken an error as a line-to-line scatter divided by the square 
root of the number of lines (e.g. an error in the mean of the 
abundances from individual lines). When only one line is 
available for the abundance estimate, then an error is taken 
to be equal to 0.25, which is the line-to-line scatter in abun- 
dances derived from individual Fe I lines. 

However, uncertainty from systematic errors may be 
much larger than the random errors from the line-to-line 
scatter. Systematic errors due to the adopted atmospheric 
parameters are checked by changing the parameters by 
ATcft = ±400 K, A log g — +0.5 dex (a model atmosphere 
with log<? = 0.0 dex is not available) and A£ = ±2.0 km 
s _1 . The results of this exercise and obtained [X/Fe] values 
are summarized in Table [3] Typically, changes in [X/Fe] val- 
ues are < 0.3 dex for the change in the T e ff and < 0.1 dex 
for the change in the logp, except for [N/Fe], [O/Fe] and 
[Fel/H]. The [X/Fe] values change less than 0.25 dex for the 
change in the £, by 2.0 km s _1 . 

Figure HJ] shows the abundance ratios ([X/Fe]) as a func- 
tion of the atomic number. 

5.3.1 CNO and Na 

The abundance pattern of CRL2688 shows enhancement 
of C, N, O relative to Fe. The obtained abundance ratios, 
[C/Fe]= 0.63 ± 0.08 and [Q/Fe] = 0.53 ± 0.25, ar e in go od 
agreement with those derived by IKIochkova et al.1 l|2000h . 

The carbon abundance is estimated from the five lines 
with the line-to-line scatter of 0.18 dex. The [C/Fe] does not 
significantly change by the change in T e ff as shown in Table 

El 

The oxygen abundance is derived from the O I line at 
6158.2 A. The abundance may be slightly overestimated 
since this line is blended with a weak Fel line at 6157.7 
A. By synthesizing a spectrum around this line for this star, 
the contribution from the Fel line is estimated to be ~ 4 
%, which is much smaller than the other errors in [O/Fe]. 
The derived abundance is also sensitive to the adopted tem- 
perature (Table |3J. Furthermore, negative non-LTE correc- 
tion of 0.1 — 0.2 dex is expected for stars with a similar 
spectral type l|Takeda fc T akada-Hidai 1998). Therefore, an 
overabundance of [O/Fe] for this star is not yet clear in the 
present study. 

The [ N/Fe] ratio of l.()2±0.2 5 is significantly lower than 
that from IKIochkova et al.l (|200Ch . who derived [N/Fe]= 2.0. 
Since as indicated in Table [3] the derived [N/Fe] abun- 
dance largely depends on the choice of parameters in the 
model atmosphere, es pecially in T e ff, the difference from 
IKIochkova et al.1 (|2000| ) may be explained by systematic er- 
rors due to the atmospheric parameters as well as by dif- 
ference in the measured EWs. Moreover, the N abundance 
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from the Nl 7468.3 A line may suffer from non-LT E ef- 



Takcda & Takada-Hidai 



Takcda & Takada-Hidai 



fects as suggested by 
non-LTE analysis by 
gests that non-LTE corrections for the abundances derived 
from Nl lines around ~ 8700 A for F-type supergiant stars 
are ~ —0.4 dex. A similar non-LTE correction is expected 
for the Nl line used in this work. Taking into account a 
possible non-LTE correction of —0.4 dex, a modest [N/Fe] 
remains, although the non-LTE calculation is needed for a 
more reliable estimate. 

The Na abundance is estimate d from two lines (568 2.63 
A and 5688.21 A). For these lines. iTakeda et all (|2003h re- 
ported that non-LTE correction to the abundances of their 
sample stars is less than 0.2 dex, which is not significant 
within the precision of the present analysis. The derived 
abundance ratio is [Na/Fe]= 0.70 ± 0.14, which is in agree- 
ment with that derived by iKlochkova et alj (|200(t ). This 
value is higher than that observed in most of the previously 
studied disk stars w ith similar metallic ity, that generally 
show the solar ratio l|Takeda et al J 120031 ) . 



o 



1 CRL 2688 

J (Up.lim) 

» IRAS 08143-4406 



20 40 60 

Atomic Number 



5.3.2 Neutron- capture elements 

The Y abundance is estimated to be [Y/Fe] = 1.07 ± 0.37 
from two lines, where the large error comes from the line- 
to-line difference in derived a bundances. This value is in 
agreement with that derived by IKlochkova et alj (|2000T i and 
significantly higher than those observed in typica l disk stars 
with similar metallicity (jEdvardsson et al"1 ri993). Such en- 
hanced [Y/Fe], h owever, is not unusual for other carbon-r ich 
post-AGB stars (|Revniers et al.ll2004l ; iReddv et al.lll997T l. 

On the other hand, the abundance of the heavier s- 
process element Ba is not as highly enhanced as that of Y. 
We note that the Ba abundance was estimated from just 
one of the three lines detected in our spectrum because the 
other two lines seem to be affected by the emission in the 
violet wing (Figure [2}. The Ba line at 5853.67 A has an 
approximately symmetric profile, which we think it is more 
reliable. For the group of the heavy s-process elements, we 
could only obtain the upper limits of the abundances. En- 
hancements similar to [Y/Fe] are certainly not observed and 
remain comparable to or lower than the [Ba/Fe] ratio. 

The ratio of heavy to light s-process elements is fre- 
quently used as a measure of the neutron exposure respon- 
sible for processing these elements. The observed [Ba/Y] = 
— 1.0 ± 0.4 is actually very low. This result has al s o bee n 
reported by the previous study of IKlochkova et all (2000). 
In Figure [7] we compared the abundance pattern of the s- 
process elements of CRL 268 8 in this work w i th tha t of other 
post-AGB stars analyzed bv lRevniers et al.l(|2004T l. Here we 
use Ba and Y as a tracer for heavy-s and light-s process 
el ements ([hs/ls1~[Ba/Y 1 and [s/Fe]~([Y/Fe] + [Ba/Fe])/2). 
InlRe^ niers et all (|2004T ). the heavy-s process abundance is 
calculated as the mean of the Ba, La, Nd and Sm abun- 
dances, while the light-s process abundance is the mean of 
the abundances derived from Y and Zr. Figure [7] shows that 
CRL 2688 has the lowest [hs/ls] compared to the other ob- 
jects. However, CRL 2688 appears to follow the same linear 
relation in [hs/ls]- [s/Fe]. 



Figure 6. Abundance pattern of CRL2688. Gray crosses and 
dotted lines show the abundance patt ern of another carbon -rich 
post AGB star IRAS 08143-4406 from lRevniers et alj 1120041 1. 



5.3.3 Alpha and Fe-peak elements 

The a and Fe-peak elements are normal compared to the val- 
ues observed in the Galactic disk stars with similar metallic- 
ity. In particular, the non-refractory element Zn abundance 
is estimated to be [Zn/Fe]= 0.04 ± 0.25 dex, which is the 
standard value observed in disk stars with similar metal- 
licity, suggesting that depletion of refractory elements due 
to condensation onto dust grains does not significantly con- 
tribute to the observed abundance pattern. Another non- 
refractory element, sulphur, shows an enhanced abundance 
ratio of [S/Fe] = 0.65±0.06 cons i stent with the previous anal- 
ysis made by IKlochkova et alj (|2000T l . This value is higher 
than typical values reported for disk dwarf/giant stars 
with s imilar metallicity studied bv lTakeda fc Takada-Hidail 
|201ll ), that generally show the solar ratio. Since the sulphur 
abundance is measured from relatively weak lines, a higher 
signal-to-noise ratio spectroscopy is needed to confirm the 
overabundance of [S/Fe]. 



6 DISCUSSION 

6.1 Circumstellar envelope 

The radial velocities of the molecular lines provide an es- 
timation of the velocity of the circumstellar outflow where 
these lines are formed. The emission components of the C2 
Phillips bands have a mean radial velocity of —38.1 km s _1 
with respect to the photospheric atomic absorption lines, 
while the emission components of the Ha, Nal (el) and 
Kl lines have a similar value of —41 km s _1 . For the ab- 
sorption components of the C2 Phillips bands, this value is 
—53.0 km s _1 , which is in good agreement with the value of 
—51.8 km s _1 derived from the CN red system. Existence 
of several velocity components of circumstellar outflows in a 
range 19-100 km s" 1 have been identified based on CO and 
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Table 3. Abundance results 



Elem 


[X/H] 


[X/Fe] 


0"linei 


N 






log g 








(dex) 


(dex) 


(dex) 




+400 K - 


-400 K 


+0.5 dex 


+2.0 km s~ x 


-2.0 km s" 1 


C I 


0.31 


0.63 


0.08 


5 


-0.13 


0.20 


0.07 


-0.01 


-0.00 


N I 


0.70 


1.02 


— 


1 


-0.41 


0.52 


0.24 


-0.09 


0.16 


O I 


0.21 


0.53 


— 


1 


-0.42 


0.56 


0.26 


-0.06 


0.09 


Na I 


0.37 


0.69 


0.14 


2 


-0.08 


0.08 


0.00 


0.00 


-0.01 


Si II 


0.11 


0.45 


0.07 


2 


-0.30 


0.33 


0.05 


-0.02 


-0.00 


S I 


0.33 


0.65 


0.06 


3 


-0.11 


0.13 


0.02 


0.02 


-0.05 


Ca I 


-0.20 


0.12 


0.14 




-0.05 


0.05 


-0.00 


-0.00 


0.00 


Sc II 


-0.01 


0.32 


0.07 


1 


0.08 


-0.07 


0.01 


0.00 


0.01 


Ti II 


-0.22 


0.11 


0.00 


2 


0.08 


-0.07 


0.01 


0.01 


0.00 


Cr II 


-0.23 


0.10 


0.07 


8 


-0.02 


0.02 


0.01 


0.03 


-0.08 


Mn I 


-0.14 


0.18 




1 


0.02 


-0.01 


0.00 


0.01 


-0.03 


Fc I 


-0.32 


0.00 


0.07 


13 


0.41 


-0.44 


-0.19 


-0.04 


0.09 


Fe II 


-0.34 


0.00 


0.04 


5 


0.25 


-0.18 


0.05 


-0.10 


0.24 


Ni I 


0.04 


0.36 




1 


-0.02 


0.03 


0.00 


0.00 


-0.00 


Zn I 


-0.29 


0.04 




1 


0.00 


0.00 


0.01 


0.03 


-0.07 


Y II 


0.73 


1.07 


0.37 


2 


0.08 


-0.08 


0.00 


0.05 


-0.10 


Ba II 


-0.26 


0.07 




1 


0.13 


-0.11 


0.01 


-0.07 


0.21 



a When only one line is available for the abundance estimate, then an error is taken to be equal 
to 0.25, which is the line-to-line scatter in abundances derived from individual Fe I lines. 
" Effect of errors in atmospheric parameters on [X/Fe] ratios. 

these molecules are confined to a thin circumstellar shell 
whose inner and outer radius are determined by the 
penetration/shields of photons from the radiation fields 
l|Bakker et al.lll997T ). The presence of both absorption and 
emission components in C2 Phillips bands suggests complex 
geometric structure of this circumstellar shell in which C2 
and CN can exist. Furthermore, the radial velocity differ- 
ences between these components indicate that the regions 
where the absorption components of the C2 and CN bands 
arise are expanding faster than those for the emission com- 
ponents. 



6.2 Abundance pattern 

The enhanced C, N and Na abundances seen in CRL 2688 
are indicative of dredge-up and mixing of H and He burning 
products during the AGB phase of evolution. In this process, 
freshly synthesized s-process elements are also expected to 
be dredged-up, giving rise to overabundance in s-process ele- 
ments in the photosphere. Indeed a class of carbon-enriched 
post-A GB stars are known to be enhanced with s-pro cess el- 
ements (|Hrivnak fc Reddvl2003l ; lRevniers et al.ll2004 ). How- 
ever, not all of them show this enhancement and some stars 
show depression in the s-process elements instead. The origin 
of this diversity i n s-process abundanc e patterns is not well 
understood (e.g. iKappeler et al.ll201ll ). Theoretical models 
for the generation of s-process elements during the AGB evo- 
lution depends strongly on the parametrization of unknown 
factors su ch as a 13 C pocke t , which may contr ol neutron 
exposure (Bus so et al. I ll999h . IBusso et all (| 19991 ) also sug- 
gested that the s-process efficiency depends on progenitor 

met allicity. 

iRevniers et alj ([2004 ) found a strong correlation of 
[hs/ls] index with the overall s-process elemental abundance 
([s/Fe]) for their sample of carbon-rich post-AGB stars. 
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Figure 7. Heavy s-process element to light s-process element 
ratios plotted against mean of s-process to Fe ratios. In this work 
(filled circle), [s/Fe] is taken as a mean of [Y/Fe] and [Ba/Fe], 
while [Ba/Y ] ratio is used as a p roxy of [hs/ls]. For the data 
obtained by IRevniers et al.l l|2004li (gray crosses) , the definition 
from the original paper is used for both quantities. 



13 CO observations at r adio wavelengths (I Young et al.lll99^ ; 
lYamamura et al1ll996l ). 

In the circumstellar environment, these molecules are 
thought to form from more complex molecules like C2H2 
and HCN that are photodissociated by interstellar and/or 
stellar radiation fields. Therefore, it is suggested that 
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They have suggested that the observed variation in the 
s-process elemental abundance patterns can partly be ex- 
plained by AGB models with varying strength of the 13 C 
pocket. The observed low [Ba/Y], treated as a proxy of 
[hs/ls], appears to follow the suggested [hs/ls]-[s/Fe] relation 
but with the lowest estimate o f [hs/ls], which is low er than 
typical theoretical predictions (|Kappeler et al-lboill ). Alter- 
natively, a neutron source of 22 Ne(a,n) 25 Mg, which would 
operate in rather massive AGB stars, might be responsible 
for the observed abundance pattern. If this is the case, the 
progenitor may be a star with more than 3 or 4 solar masses. 
However this would be in contradiction with the observed 
C-rich nature of the circumstellar envelope of CRL 2688. 
Note that the abundance result suggests that C/0<1 from 
the photospheric absorption lines. However, the oxygen lines 
may be influenced by non-LTE effects. If we apply non-LTE 
correction to oxygen abundance the C/O ratio my be 1.0 or 
more. 

Since only Y and Ba have reliably been estimated in 
the present work, the conclusion about the very low [hs/ls] 
remains to be confirmed. Further accurate abundance study 
with much higher resolution and high signal to ratio spec- 
trum may enable us to understand nucleosynthesis in the 
progenitor star. 



7 CONCLUSION 

We have presented the analysis of a high-resolution opti- 
cal spectrum of CRL 2688. Based on the one-dimensional 
LTE abundance analysis, we have confirmed that the cen- 
tral star has a spectral type of a F-type supergiant with 
[Fe/H]= -0.3±0.1. Molecular absorptions from C 2 Phillips, 
C2 Swan bands and CN red system, that are presumably of 
circumstellar origin have been identified, some of which ap- 
pear in emission. In particular, C2 Phillips band show pres- 
ence of emission components with the helio-centric radial 
velocity offsetted from the photospheric absorption lines by 
—38 km s , which is consistent with emission components 
seen also in the Ha line, Nal and Kl resonance lines (—41 
km s _1 ). The absorption components of these lines have a 
mean radial velocity of —52 km s _1 with respect to the pho- 
tospheric absorption lines, and are more highly blueshifted 
than the corresponding emission components. 

The abundance pattern of CRL2688 is characterized by 
the enhancement of C, N, Na and Y. However, unlike other 
known C-rich post AGB stars, this object shows a very low 
[Ba/Y] ratio. Since only abundances of Ba and Y have been 
reliably measured as a tracer of heavy and light s-process 
elements, respectively, in the present work, the conclusion 
about the peculiar [hs/ls] ratio observed in CRL 2688 re- 
mains to be confirmed with measurements of more heavy 
and light s-process elements. 

Higher resolution and higher signal-to-noise ratio spec- 
troscopy is required for more comprehensive understanding 
of the origin of the each spectral feature and the detailed 
abundance pattern. 
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